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Abstract We study two common types of time-noncontinuous updates 
for one dimensional stochastic cellular automata with arbitrary nearest neigh- 
bor interactions and arbitrary open boundary conditions. We first construct 
the stationary states using the matrix product formalism. This construction 
then allows to prove a general connection between the stationary states which 
are produced by the two different types of updates. 

Using this connection, we derive explicit relations between the densities and 
correlation functions for these different stationary states. 



1 email: nr@thp.uni-koeln.de 

2 email: schreckOmath . uni-duisburg . de 



In this letter we study one dimensional stochastic cellular automata with 
different types of update. We consider a chain of sites i — 1, ...,L (L even) 
and assign a discrete variable Tj = 1, m to each site i. Tj can for example 
represent m different sorts of particles. We then prescribe some local rule 
T describing the interaction between a pair of nearest neighbor variables 
(Tj,Tj + i) and, as usually done, choose this rule to be independent of the lo- 
cation of the sites. One can further allow for some rules TZ and C, acting on 
the sites at the two ends of the chains. Such boundary conditions could for 
example mimic the injection and removal of particles with certain rates. This 
type of boundaries can lead to a rich dynamical behaviour of the model. One 
prominent example is the asymmetric exclusion model (see U, || and ref. 
therein) which exhibits boundary induced phase transitions where the bulk 
density of particles depends in a subtle way on the boundary conditions. 
There are now several ways to implement our prescribed update rules. These 
implementations are an essential part of the definition of a model since they 
generally produce different transients, stationary states and averages of ob- 
servables. 

For analytic calculations one usually uses the random sequential update. One 
picks at random site i and applies the update rule. This update generally pro- 
duces the weakest correlations between the sites compared to other update 
procedures. The master equation which describes the temporal evolution 
of the weights of the configurations {t\,T2, ■■■,tl} can then be written as a 
Schrodinger-like equation with a (non-hermetian) hamiltonian consisting of 
a sum of local hamiltonians. For this update it was recently shown, that the 
stationary state can be always written as the scalar product of the product 
of some (generally) noncommuting matrices 0. 

There are three other basic ways of implementing the update rules which are 
much more common for practical purposes like computer simulations, since 
the random sequential update requires an additional random number for each 
local update. These three types are called sequential, sublattice parallel and 
parallel update. In the parallel update one applies the update rules simulta- 
neously to all sites. This case will not be discussed in the following. Note 
that the parallel update usually produces the strongest correlations between 
the sites and is often used in the context of traffic flow models [|]. 
In the sequential update we start at the right end of the chain and apply 
our boundary condition 71 to site i = L. We then update the pair (rx_i, tl). 
After that we proceed with the pair (tl_2, Tl-i) and so forth until we reach 
the left end of the chain where we apply £ at site i = 1. This procedure can 
of course also be done from the left to the right. Let us denote these two 
updates with and 

The sublattice parallel update 21 1 is defined as follows: In the first time step 



2 



one applies C, 1Z and updates all pairs (rj, Tj+i) with an even i. In the second 
time step one updates all pairs (7$, r i+1 ) with an odd z. 
In the following, we show that in the stationary state averages like densities 
and correlation functions of the sequential and sublattice parallel update are 
essentially the same. This is a nontrivial result since the microscopic config- 
urations of the system are generally very different after applying a sequential 
or a sublattice parallel update. This result should be of practical use for 
people doing computer simulations of cellular automata or reaction-diffusion 
processes. Note that the sequential update can transport a particle many 
sites during one update of the whole chain while T» can move a particle only 
two sites at maximum. 

We start with the construction of the stationary state of an arbitrary model 
with sequential update T<_ from the right to the left end. We will first show 
that it is always possible to write the stationary probability distribution 
-Po( T "i) 7 "2) ■••jTl) of the configurations {r^} as 

Po(t 1 ,t 2 ,...,t l )^ = (W\A Ti A T2 ...A Tl \V) (1) 

where the A T . are matrices (representating single site states), the (W\, \ V) 
vectors (reflecting the influence of the boundary conditions) acting in some 
auxiliary space. One can formally define a column vector A = (A\, ...,A m ) 
and rewrite ([!]) as 

|P >- = (W\A 9L \V) . (2) 

The existence of a nontrivial representation of (^) can then be shown simply 
by defining the AiS, (W\ and | V) exactly as in ref. 0. We use here a slightly 
different notation and define the matrix A T (r = 1, 2, m) by 

A = \t)(V\ + J2\tt2)(t 2 \+ £ \tt 2 t 3 )(t 2 t 3 \ + ... . (3) 

{T2} {T2TS,} 

The vacuum vector is denoted by \V). This definition implies 

A T1 A T2 ...A TL \V) = \t 1 t 2 ...t l ) . (4) 

One can then choose (W\ such that (H) holds by definition. 
Since |po)<- is defined as the stationary state one has T^\p )^ = |p )<— • The 
simplest mechanism to achieve this is to imagine that TZ produces a 'defect' 
A which is transported through the sequential action of T to the left end to 
the chain where it disappears by applying C. This mechanism reads then: 

T [A <g> A] = A g> A , (5) 

n[A]\v) = A\v) , (6) 

(W\C[A} = (W\A. (7) 
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By defining 



with 



L-1 



A® A® \V) = J] TiRA 9 \V) 



8=1 



R = 1®...®1®K, (9) 
Ti = 10 1...®%® (10) 

it can be shown analogously to |j| that in fact the defined A fulfills the 
algebra©-®. 

This means that in order to solve such a model exactly, the problem of solving 

the master equation has been reformulated in terms of finding representations 

of the algebra. Up to now, explicit representations of this algebra are only 

known for the asymmetric exclusion model ( ||,|J,0 )• 

We are now able to write down the stationary states for the updates and 

T f 

\PoU = (W\A® L \V) , (11) 

and 

|P )ll = (W\ A® A® A® A® ... <g> A® A\V) . (12) 

By applying Tl* on ([□]), T\\ on ([T2J) and using the algebra (§j)-(0) one can 
confirm that these states are in fact the desired stationary states. One still 
has to confirm that these expressions are not equal to zero. To see this we 
define 

m m 

C = Y,A ; <? = EA • ( 13 ) 

1=1 1=1 

The matrices T,7Z,£ have the property that their columns add up to one since 
the update of a local state has to result in some other local state. Adding 
up each equation (|§)-(0) one obtains 

[C,C] = 0, (14) 
C\V) = C\V), (15) 
(W\C = (W\C. (16) 

This implies immediately that the norm (W r |C i |V r ) of |p )<— is equal to the 
norms of |po)— >• anc i \Po)\\ which shows that the latter states are nontrivial. 
One can now compare the densities of the stationary states produced by the 
different types of update. Using the formulas <g), QTT]), (0) and (^-(03) 
one gets for the densities Pi{x), which are defined as the probability to find 
a particle of type r — 1, m at location i = 1, L in the stationary state: 

P A^ = l P ih : lodd ■ (it) 

w ll I p T [i)^ : zeven 
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Note that reversing the order of the the two time steps in the sublattice 
parallel update would result in reverse arrows in (|17D . 

In order to get similar expressions for the correlation functions we now make 
use of the identity 

\v)(v\+ J2 E I^^-^iXti,^,...,^! = i, (is) 

1=1,2,... {-n,T2, — ,Tl} 

which expresses the completeness of our Fock space. 
This can be rewritten as 

E C l \V)(V\{C l ) ] = 1 . (19) 

{=0,1,2,... 



One can now multiply both sides of fliTf ) with C and make use of ([TJ]) and 
(p~5f ) which leads to 



C = C. (20) 

This result can actually help to construct explicit representations of the alge- 
bra since it reduces the number of matrices 2m which have to be representated 
by one. It also reduces the m 2 equations for the bulk algebra by one and the 
2m equations for the boundary algebra by two, which can be of significant 
help. 

Note that in the case of random sequential update the same considerations 
would lead to the result (7 = 0, which is fulfilled by the matrix product so- 
lution of Derrida et al of the asymmetric exclusion model ||. 
One further obtains a relationship between the correlation functions in the 
stationary states which reads for the two-point correlation function (gener- 
alizations to n— point correlation functions are obvious): 

(r.r,),, ={ < T ' T '->- : ^:° dd (21) 
x W 1 {TiTj}^_ : i, j even 

The correlations between odd-even and even-odd sites cannot be expressed 
so easily. Since the matrix C acts like a transfer matrix between the points 
of the correlation functions, the relevant length scales are determined by C 
and will therefore be equal for all three updates. 

Since the time evolution of a system can be written as \p t ) = T'|p i=0 ) it is 
useful to ask for the general eigenvectors \p E ) with eigenvalue E of T. It turns 
out that these can also be constructed (|[L0||,p|) by writing E = e l eie r and 
multiplying the r.h.s of (|8p with eL -i e ■ Details will be given in [jSJ. It turns 
also out that by setting e = 1 and e\ = j- one has the freedom to substitute 
I) with C = ±C. 
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In this letter we have shown that the matrix product formalism can be used 
to prove the physical equivalence of the sequential and the sublattice parallel 
update. One gets explicit nontrivial connections between densities and cor- 
relation functions in the stationary state. 

One can make similar considerations concerning the relation between the 
stationary states of the random sequential update and the updates discussed 
in the present work. This work is in progress f|. 
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